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ABSTRACT 

We have used the Far Ultraviolet Spectroscopic Explorer to search for O VI A1031.926 emission at 
four focations in the star burst superwind of M82. No O VI emission was detected at any of the four 
pointings, with upper hmits less than or equal to the 0.3—2 keV X-ray flux. These observations limit 
the energy lost through radiative cooling of coronal phase (T ^ 10^'^ K) gas to roughly the same 
magnitude as that lost in the hot phase through X-ray emission, which has been shown to be small. 
The wind material retains most of its energy and should be able to escape from the gravitational 
potential of M82, enriching the intergalactic medium with energy and metals. The lack of coronal 
gas in the wind and observations of spatially correlated X-ray and Ha emission are consistent with a 
scenario in which the hot wind material over-runs cold clouds in the halo, or one where the Ha and 
X-ray emission arise at the interface between the hot wind and a cool shell of swept-up ISM, as long 
as the shock velocity is < 150 kms~^. The observed limits on the O VI/Ha and C III/Ha flux ratios 
rule out shock heating as the source of the T — 10'* K gas unless the shock velocity is < 90 kms^*. 
Subject headings: Galaxies: individual (M82) — galaxies: starburst — galaxies: halos — galaxies: 
ISM 



1. INTRODUCTION 

A common product of starbursts in galaxies is a 
galactic-scale outflow of gas and energy from the star- 
burst region, called a starburst superwind (Heckman, 
Armus, & Miley 1990; Strickland & Stevens 2000). A 
superwind begins as an expanding bubble of hot gas pow- 
ered by the combined energy of stellar winds and super- 
novae. This bubble may break out of the disk, resulting 
in a wind of hot (T ^ 10^ K), low-de nsity gas flowing 
into the halo fChevahe r fc Cleg-e^ '1985l. As this hot gas 
collides with cold ambient material in the halo, interme- 
diate temperature gas (T ~ 10^ — 10^ K) is formed at the 
interface. This gas emits X-rays, and if it cools through 
T - lO^-s K it win produce O VI AA1031.926, 1037.617 
line emission. 

The question of whether superwinds are able to over- 
come the gravitational potential of their host galaxies 
is critical for determining their full impact on the evo- 
lution of galaxies and the intergalactic medium (IGM). 
If the wind does escape, the wind material will be de- 
posited into the IGM, enriching it with energy and metals 
{e.g., Aguirre et al. 2001), and it may produce the ob- 
served mass-metallicity relationship in galaxies if winds 
are preferentially able to escap e from low mass systems 
ijHeckman et al.l l2000t lMartinl[T999(l . A comparison of 
wind velocities with escape velocities suggests that winds 
can readily escape fro m dwarf galaxies, w hile winds of 
larger galaxies cannot l|Heckman et al.ll20b0^') . However, 
wind velocities have only been measured relatively close 
to the starburst region, so it is still unclear whether en- 
ergy loss further downstream keeps the material bound 
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to the galaxy. Observations and models imply that ra- 
diative cooling through X-ray emission from hot gas is 
not effective at removing kinetic energy from superwinds 
ijStrickland fc S tevens 2000). However, coronal tempera- 
ture gas at T ~ lO^ '^ K cools very quickly and efficiently 
({Sutherland fc DopitallT99 31 . so energy lost through this 
phase may be enough to prevent winds from escaping. 

The O VI doublet at 1031.926, 1037.61 7 A is the most 
impo rtant coolant for coronal gas (Edg ar fc Chevalieil 
Il986j) . making O VI emission an excellent measure of 
the energy lost in the coronal phase. This wavelength 
range is inaccessible from ground-based telescopes, and 
opportunities to observe it from space have been limited. 
This has changed with the launch of Far Ultraviolet Spec- 
troscopic Explorer (FUSE; Moos et al. 2000). We have 
used FU SE to search for O VI emission in the halo of 
M82, the prototypical starburst superwind galaxy. 

2. observations and data reduction 

We obtained FU SE spectra at four positions in the 
halo of M82 on 2002 February 3. The coordinates and 
exposure times are listed in Table 1, and the positions of 
the pointings can be seen in Figure 1. The large (LWRS: 
30"x30") apertures were used. FUSE consists of four 
separate telescopes, two with LiF coated optics for op- 
timal sensitivity to longer wavelengths (1000 to 1180 A) 
and two with SiC coated optics for optimal sensitivity to 
shorter wavelengths (900 to 1000 A). These are usually 
referred to as the LiFl, LiF2, SiCl, and SiC2 channels. 
Of these, the LiF channels have the largest effective area 
at 1032 A. 

The raw photon event lists were screened to remove 
data taken while the satellite was in the South Atlantic 
Anomaly, and to remove events with pulse height ampli- 
tudes less than 4 and greater than 15, because these are 
predominantly background events. This photon list was 
then transformed to a 2-dimensional image from which 
we extracted the spectrum, using bright airglow lines to 
determine its location on the detector. The extracted 
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Fig. I.— Left: An Ha image of M82 ICheng. et al.ll996l) . showing the locations of the FUSE pointings. Right: The Cftandra 0.3 - 2.0 keV 
image of M82 (Strickland et al. 2003). The ellipse shows the approximate location of the optical disk of M82. 



Table 1. FUSE Observations 



Pointing 


R.A. 


Docl. 


Exposure Time 




(J2000.0) 


{J2000.0) 


{^) 


A 


9 55 57.1 


69 39 49.5 


8056 


B 


9 55 60.0 


69 39 17.1 


8768 


C 


9 56 03.1 


69 38 48.0 


11529 


D 


9 55 52.3 


69 42 54.5 


10011 



spectra from the individual exposures were co-added to 
produce the final spectrum for each pointing. The wave- 
length solution from the CALFUSE pipeline (v2.0.5) was 
applied to the data, but we kept the data in units of 
counts. The individual channels were not co-added. 

3. LIMITS ON O VI EMISSION 

No O VI A1031.926 or A1037.617 emission was evident 
in any of the spectra. We also searched for C III A977.020 
emission in the SiC channels, finding none. We base 
our O VI limits on the 1031.926 A line because it is the 
strongest of the O VI doublet. The background flux level 
was determined by taking the mean of the flux measured 
in ten 0.34 A bins near 1032 A. The noise in the back- 
ground is the primary source of uncertainty in the O VI 
measurement, so the square root of the background flux 
corresponds to the Icr uncertainty, from which we derive 
3(7 O VI detection limits of the FUSE spectra. 

The choice of 0.34 A bins for computing the sensitiv- 
ity corresponds to the filled aperture resolution for the 
LWRS aperture, and also matches the expected veloc- 
ity range of O VI emission based on absorp tion obser- 
vations of NGC 1705 ffleckman et alJl200l() . However, 
IShopbell fc Bland-Hawthornl l|1998(r found that the Ha 
emission is spread over 300 kms-i in the M82 outflow, 
which would correspond to an emission FWHM of 1 A at 
1032 A. Using this bin size would increase the observed 
O VI upper limits by a factor of 1.7. We adopt the for- 
mer value (FWHM= 100 knis^^) since it is based on an 



O VI measurement, but the difference is not enough to 
change our conclusions. 

These limits were converted from counts to flux us- 
ing the effective area of the LiFl channel at 1032 A and 
the exposure time. The observed limit was corrected 
for G alactic extinction of Ay =0.10 (Burstei n fc Heile3 
Il984|) and extinction intrinsic to the M82 superwind 
Ay = 0.85 iHec kman et al. 192Q|). We also estimated 
the extinction by modeling the photoel ectric absorption 
in Ch andra 0.3 — 2.0 keV X-ray data ijStrickland et alJ 
12003(1 ■ and found agreement with the extinction de- 
rived from the Balmer decrement by IHeckman et alJ 
(1990). The X-ray opacity is dependent upon dust grains 
and the gas-phase metal column, so it is directly re- 
lated to the ISM constituents responsible for the far- 
ultraviolet extinction. We corrected the internal extinc- 
tion three different ways: using the Galactic extinction 
law (CardeUi. Clayton, & Mathis 1989), the LMC ex- 
tinction law (Howarth 1983 ), and the starburst extinc - 
tion law ( Calzetti. Kinnev. fc Storchi-BergmannI [1994(1 . 
The LMC law produces more far-ultraviolet extinction 
than the Galactic law, while the starburst law produces 
much less. We adopt the Galactic extinction law as the 
intermediate case, although using any of the three laws 
would not alter our main conclusions. The extinction- 
corrected limits are listed in Table 2. Figure 2 shows the 
observed spectrum of position A (see Figure 1), along 
with a model of the expected O VI emission (if it were 
equal to the X-ray flux). 

Position D is on the North side of the M82 outflow, and 
since the sight line passes through the disk of M82 it is 
likely subject to higher extinction. The optical spectrum 
of Heckman et al. (1990) did not reach this far into the 
halo, so we have no measure of the extinction in this 
region. The most reliable limits come from positions A, 
B, and C, which are on the southwest (near) side of the 
outflow and do not pass through the M82 disk. 

The velocity of M82 [vr = 203 kms^^, de Vaucouleurs 
et al. 1991) shifts the O VI emission to 1032.624 A. The 
only Galactic absorption nearby is the R(4) 6 — line 
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Table 2. Measured Properties^ 



Pointing 


Fio32 Limit'^ 
(erg cm~^ 


^coronal Limit*^ 

(erg cm~^ ^~^) 


^ coronal/^ X — ray 

(erg cm— ^ ^~^) 


Fc III Limitb 




Fc iii/Fhc," 


A 


< 7.6 X 10^14 


< 3.8 X 10-13 


< 0.40 


< 4.4 X 10-13 


< 0.01 


< 0.05 


B 


< 7.6 X 10-14 


< 3.8 X 10-13 


< 0.54 


< 4.4 X 10-13 


< 0.03 


< 0.18 


C 


< 6.2 X IQ-i* 


< 3.1 X 10-13 


< 2.58 


< 3.2 X 10-13 


< 0.06 


< 0.29 


D 


< 6.3 X 10-" 


< 3.2 X 10-13 


< 5.08 


< 3.5 X 10-13 


< 0.02 


< 0.09 



"Upper limits are 3(7. 

''Corrected for extinction, assuming a Galactic extinction of A^=0.10 magnitudes, and an M82 extinction of Av=0.85 (Heckman et al. 
1990), and using the Galactic extinction law (Cardelli et al. 1989), which gives A1032 = 4.84Ay (assuming ijy=3.1). The LMC extinction 
law (Howarth 1983) gives A1032 = 5.55j4v, which would result in upper limits ~ 1.8x greater, and the starburst extinction law (Calzetti 
et al. 1994) gives j4io32 = 2.35^4^ (normalized to match the Galactic law at 5500A), which would result in upper limits ~ 7. Ox lower. 

'^Total energy radiated by the coronal gas, corrected for extinction. 

''The absorption corrected X-ray flux (0.3-2.0 keV) within the 30" X 30" LWRS aperture was measured on Chandra ACIS-S spectra 
(Strickland et al. 2003). 

"^The extinction corrected Ha flux was measured in the 30" X 30" LWRS aperture using the image shown in Figure 1. 



of H2 at 1032.351 A. An examination of sight lines near 
M82 (~ 7° away) observed by FUSE suggests that H2 
in the J — 4 rotational level is present in this direction. 
This line is ~ 80 kms""'^ from the expected center of 
the O VI line, so the contamination is expected to be 
minimal even if the lines are broad. The M82 C III line 
falls at 977.681 A, which is in a region free of Galactic 
absorption, but is 85 kms ~^ from an atmosph eric O I 
emission line at 977.959 A ijFeldman et al.ll200lD . 

In the optically thin case, the energy radiated through 
the two O VI lines is 150% of the 1032 A line strenth. 
The two O VI lines together are responsible for ^ 30% 
of the cooling in coronal gas l|Heckman et al.l [2001). The 
limit on the total energy radiated by coronal gas Fcoronai 
is thus a factor of 5 higher than the limit on 1032 A 
flux (see Table 2). These limits are comparable to the 
observed X-ray flux at each position, ranging from 40% 
to 260 % in positions A, B, and C. Strickland fc Stevens! 
l)2000tl found that the energy lost through X-ray emission 
in M82 is < 10% of the energy input from the starburst. 
If the lowest upper limit on O VI is valid everywhere in 
the wind, then the amount of energy radiated away by 
the coronal gas is even smaller than that lost through 
X-rays. 

4. DISCUSSION 

4.1. Radiative Cooling in the Superwind 

The non-detection of O VI emission has profound im- 
plications for the evolution of the starburst superwind 
in M82. The lack of O VI emission indicates that ra- 
diative cooling from coronal gas docs not effectively re- 
move energy from the wind. Energy loss through other 
phases has been deter mined to be insignificant as well 
ijStrickland fc Stevens! 120001 . so the wind must retain 
most of its kinetic/thermal energy as it travels through 
the halo. The wind velocity comfortably exceeds the es- 
cape velocity, so the wind should easily escape from the 
galaxjr^ 

Ot te et al] l)2003f ) detected O VI emission in the halo 
of NGC 4631 at the locations of strong X-ray emis- 
sion. NGC 4631 is actively forming stars but is not a 
starburst and has weaker X-ray emission overall than 
M82. Nevertheless, the O Vl/X-ray ratios were ^ 1, 
similar to our limits in M82, suggesting that radiative 



cooling of coronal gas is ineffective in the hot halo of 
NGC 4631, a conclusion we also reach for the M82 super- 
wind. Heckman ct al. ( 2001) measured O VI absorption 
and set upper bounds on the O VI emission in the star- 
burst NGC 1705, and were able to limit the cooling rate 
to < 20% of the supernova heating rate in that galaxy. 
Since the X-ray luminosity of NGC 1705 is very low, the 
NGC 1705 superwind should also escape into the IGM. 
If these results hold for many starbursts, they could be 
an important source of energy and metals for the IGM. 

4.2. Origin of O VI, Ha, and X-rays 

IStrickland et al.l lj2002j) found significant spatial cor- 
relation between Ha and X-ray emission in the star- 
burst superwind of NGC 253, an d a similar correla- 
tion is seen in the M82 superwind ijLehnert et alJir999l: 




-1 P. . . I I I I . 

1031 1032 1033 1034 

Wavelength (A) 



Fig. 2. — The FUSE spectrum of the wavelength region around 
1032 A in position A. The observed spectrum (histogram) is binned 
to 0.037 A/pixel, and a constant background was subtracted. No 
O VI emission is visible in the spectrum. The smooth lines show 
the expected number of counts if the O VI emission were equal to 
the observed 0.3-2.0 keV X-ray flux, after accounting for extinc- 
tion using the Galactic extinction law (Cardelli et al. 1989, solid 
line), and the LMC extinction law (Howarth et al. 1983, dashed 
line). The FWHM of the emission line is assumed to be 0.34 A 
(100 km s-l). 
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iStrickland et al.1 l2003j) . This fact strongly suggests a 
physical connection between the origin of both types of 
emission. Models of super winds must account for this 
correlation, and [Strickland et al.l l)2002(l outline several 
possibilities. The FUSE pointings were chosen because 
they exhibit strong X-ray emission and Ha emission, 
suggesting the presence of warm and hot gas (see Fig- 
ure 1). 

It is unlikely that both the X-ray and Ha emission 
arise in cooling wind material. IStrickland et al.l l|2002fl 
showed that the cooling times are too long, and the 
non-detection of O VI emission indicates that significant 
amounts of gas do not cool through T = 10^'^ K. Another 
possibility is that the Ha emission arises in a cool shell of 
swept up ISM surrounding the superwind {e.g. Weaver 
et al. 1977). If the shock velocity Vs > 150 kms~^ the 
interface between the hot wind fluid a nd the cool shell 
mater ial would contain coronal gas, and IStrickland et al.l 
l)2002D showed that the luminosity of the O VI doublet 
would be ~ 100 times greater than the X-ray luminosity. 
This is clearly at odds with the observed O VI upper lim- 
its, but we cannot rule out a lower Vs- We also cannot 
rule out a hot shell of swept-up ISM which would not 
produce coronal gas due to long cooling times, but this 
scenario has difficulty reproducing the spatial correlation 
betw een the X-ray and Ha emission (Strickland et 

l2pL ,__ 

ILehnert et al.l l|1999|) discuss a model in which the wind 
runs into cool clouds in the halo, either pre-existing or 
entrained within the wind itself (see also Strickland et al. 
2002) . The X-ray emission is produced in a stand-off bow 
shock upstream from a cloud, while the Ha is produced 
by a slow shock driven into the cloud by the wind (with 
a possible contribution from photoionization). The X- 
ray and Ha emission in this scenario are physically cor- 
related, but little coronal temperature gas is produced, 
consistent with our non-detection of O VI emission. 

Comparing the O VI and C III limits with the Ha 
fiux can also constrain the conditions in the super- 
wind. If a turbulent mixing layer between hot and 
cold gas forms in the wind, the intermediate temper- 
ature gas produced must be cooler than 10^'^ K, or 
Fo vf/Fff„ would be higher than the obse rved limits 
in Table 2 l|Slavin. Shull. fc BegelmM]IT99l . Figure 3 
shows the predictions of shock models for the O VI/ Ha 
and C III/ Ha flux r atios (Dooita & Sutherland 1996^ 
iShull fc McKeel n.979|^ . compared to the observed up- 
per limits on the ratios at position A. The upper limits 
shown were calculated assuming a line width (FWHM) 
corresponding to the shock velocity, so the limits are 
higher for broader lines. The limits for all four point- 
ings are listed in Table 2, assuming the line width 
is 100 kms~^. The Fo vi/^Ha upper limit is lower 
than the prediction for shocks with Vs ^ 150 km s^^, 



and the Fc iii/^Ha ratio is lower than the predic- 
tion for shocks with Vs ^ 90 km s~^, assuming all of 
the Ha emission arises in shock-heated gas (very fast 
shocks are allowed by the Fc iii/Fhch upper limit, but 
this is ruled out by the Fn yj / F Hn upper limit). If 
the model of ILehnert et al.l l)1999.il is correct, the slow 
shocks in the cold clouds that produce Ha emission have 
Vs ^ 90 km s~^, or the ionization of the cold clouds 
is dominated by photoionization from the starburst. 
CLOUDY models ijFerland et all Il998j) show that the 
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Fig. 3. — Shock model predictions of O VI/Ha (solid line) 
and C III/Ha (dashed line) from the models DoDita &; Sutherlandj 
11993) (vs > 150 kms-l) and Shull & McKcc 11973) i^a < 
150 kms~^). The nearly diagonal lines show the observed upper 
limits at position A, the most restrictive limit of the four pointings. 
The upper limits are calculated for a line width (FWHM) corre- 
sponding to the shock velocity (assuming bulk motion, not thermal 
broadening, is responsible for the line width). The upper limits are 
well below the predictions, indicating that most of the Ha-emitting 
gas is not shock heated, unless the shocks have Vs < 90 km . 



Fc III /F Ha limit is consistent with photoionization (this 
will be discussed further in a future paper, C. Hoopes et 
al., in preparation). If the clouds are photoionized, this 
means that the energy source powering the Ha emission 
is not the wind, so even the optical emission does not 
represent energy lost from the wind fluid. This scenario 
raises the possibility that the wind blown cavity provides 
a conduit for ionizing photons from the starburst to reach 
the IGM. 

We thank Ken Scmbach for helpful discussions. We 
appreciate the helpful comments from the anonymous 
referee. This work was supported by NASA grant NAG5- 
11945. 
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